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Abstract: In this work we generalized holographic model for s-wave DC Josephson 
junction constructed in arXiv:1101.3326[hep-th] to a holographic description for p-wave 
Josephson junction. By solving numerically the coupled equations of motion of Yang-Mills 
theory for a non-Abelian SU(2) gauge fields in (3+l)-dimensional AdS spacetimes, we 
shown that DC current of the p-wave Josephson junction is proportional to the sine of the 
phase difference across the junction like the s-wave case. 
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1 Introduction 

The AdS/CFT correspondence [1] has attracted much attention over the past few years. 
The AdS/CFT correspondence relates a gravity theory in a weakly curved (d + l)-dimensional 
AdS spacetimes to the strong coupled d-dimensional field theory living at the AdS bound- 
ary. In the last few years, with the application of the AdS/CFT correspondence to the 
superconducting phase, ones find that black holes with charged scalar hair in AdS space- 
times can provide a holographically dual description of superconductivity [2-4] . For reviews 
on holographic superconductors, see [5-7]. 

Recently, a holographic model for s-wave Josephson junction with DC current con- 
structed in [8]. Solving numerically the equations of motion of Einstein-Maxwell-scalar 
model with spatial dependence \i and non-vanishing constant J, ones show that the DC 
current of s-wave Josephson junction is proportional to the sine of the phase difference 
across the junction. Meanwhile, they also studied dependence of the maximum current 
on the temperature and size of the junction, which matches precisely with the results 
for condensed matter physics. This holographic model is extended to (3+l)-dimensional 
Josephson junction in [9, 10]. Another different way to construct a holographic Josephson 
junction based on designer multigravity is studied in [11]. 

As is well known, in condensed matter systems, SNS Josephson junction can be made 
up of two p-wave superconductors. A holographic description of p-wave superconductors 
has been developed in [12-14]. Naturally, it would also be of interest to set up a holo- 
graphic model for p-wave Josephson junction. In this paper, we would like to generalize 
the work [8] to a holographic description for the p-wave Josephson junction. We will study 
a non-Abelian SU(2) Yang-Mills theory in (3+l)-dimensional AdS spacetimes and solve 
numerically the coupled, nonlinear partial differential equations of motion, and analyze the 
phase dependence of the Josephson current between two p-wave superconductors. 

The paper is organized as follows. In Sec. 2, we review non-Abelian SU (2) Yang- 
Mills theory in (3+l)-dimensional AdS spacetimes and set up a gravity dual of a (2+1)- 
dimensional p-wave Josephson junction. In Sec. 3, we show numerical results of the EOMs 
and study the characteristics of the (2+l)-dimensional holographic p-wave Josephson junc- 
tions. The last section is devoted to conclusion. 
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2 Holographic setup of p-wave Josephson junction 

Now let us concentrate on the dynamics of Yang-Mills theory for a non-Abelian SU (2) 
gauge fields in (3+l)-dimensional AdS spacetimes. First, we write the action of (3+1)- 
dimensional AdS gravity as 

S = j d 4 x^(R+^), (2.1) 

where L is the curvature radius of asymptotic AdS Spacetimes. The planar Schwarzschild 
black hole can be written as: 

ds 2 = -f(r)dt 2 + ^- + r 2 (dx 2 + dy 2 ) , f(r) = ^ - — , (2.2) 
j[r) r 

where M is the mass of the black hole. The Hawking temperature of the black hole reads 
as 

^_ 1 df 



4tt dr 



4vr L 4 /3 

r=ru 



3 M ' /3 (2.3) 



I 2 

where ru = M'iL'i is the event horizon of the black hole. In this background of gravity, 
we now consider a non-Abelian SU (2) gauge fields, with the action 



S= I d*x^[- 7 F^F^), (2.4) 



1 

? 

where F£ v = d^A^ - duA^ + e abc A b ^A c u is the field strength of the SU(2) gauge theory, 
and one form A = A^dx^ = A a ^T a dx^ is gauge field, r a are the generators of SU(2) with 
a = 1,2,3. 

Variation of the action (2.1) with respect to the gauge field A^ leads to the equations 
of motion 

D^F/iv = , (2.5) 

where the gauge covariant derivative is = V M + i [A^ , ] . 

In order to obtain a direction along which the SNS stack is arranged, we take the gauge 
field ansatz as 

A = <p(r, y)r 3 dt + w(r, y) r l dx + A r (r, y) r 3 dr + A y (r, y) r 3 dy , (2.6) 

where <fi, w, A r , and A y are real functions of r and y. Thus the junction can be along the 
y direction. 
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With the black hole background (2.2) and the above ansatz (2.6), the equations of 
SU(2) gauge theory (2.5) can be written as: 



d r w + ^J d yy w + j9 r w+ [^-A 2 r -^ I )\ij\=0, (2.7) 




= , (2.8) 
= , (2.9) 
= , (2.10) 

= 0. (2.11) 

where a prime denotes derivative with respect to r. Because Eqs. (2.7)-(2.11) are cou- 
pled nonlinear equations, one can not solve these equations analytically. However, it is 
straightforward to solve them numerically. 

In order to solve these coupled equations, first, we need to impose regularity at the 
horizon and the boundary on the radial coordinate. Near the horizon (r = r#), the field (f> 
should be regular: 

<t>{r H ) = (2.12) 
Near the AdS boundary (r — > oo), the Yang-Mills fields take the asymptotic forms 

. = + ^ + (2,3) 

* = M(l/)-^ + o(^) , (2-14) 
A r = O (J^j , (2.15) 
Ay = v(y) + ^ + o(J^ , (2.16) 

here, (i is the chemical potential and p is the charge density, v and J are the superfluid 
velocity and current, respectively [15, 16]. It is worth pointing out that current J is a 
constant and along the y direction. For u>W and are normalizable, one can impose the 
condition either w/ 1 ) or vanishes. For simplicity, we will adopt the constraint u;W = 
and interpret (O) = as the condensate. Moreover, we introduce the phase difference 
7 = — j Aydy, which is the gauge invariant and can be rewritten as 

/oo 
dy [u(y) - u(±oo)} . (2.17) 
-oo 

Furthermore, we impose the Dirichlet-like boundary condition on A r and Neumann- 
like boundary condition on w, 4> and A y at the spatial coordinate y = 0. At y = ±oo, the 
field functions are y-independent. Thus the boundary conditions of the coupled equations 
(2.7)-(2.11) are determined by J and /j,. 
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(a) 4, (b) A y 

Figure 1. The components cf> and A y of Yang-Mills fields. The parameters are set to J/T^ = 0.0094, 
A*oo = 5, L = 3, e = 0.6, a = 0.5. 



Ones need introduce the critical temperature T c of the p-wave junction, which is pro- 
portional to /x(oo) = fi(-oo) : 

T c = i-/fM > (218) 

47T I2 C 

where [i c ~ 3.65. The effective critical temperature inside the gap can be written as 

T = f^. (2.19) 

47T fl c 

In order to describe an SNS Josephen junction, we also choose the below fi(y) which is the 
same as that in [8]: 



1 - e 



2tanh(£; 



tanh 



y + k 



tanh 



(2.20) 



where the chemical potential \i is proportional to [i^ = /i(oo) = fi(-oo) at y = ±00, and L 
is the width of p-wave junction. The parameters a and e control the steepness and depth 
of p-wave junction, respectively. 



3 Numerical results 

In this section, we will solve the coupled equations (2.7)-(2.11) numerically with spectral 
methods. In order to solve EOMs numerically , we need set the change of variables z = 
1 — rn/r and y = tanh(^). First, the graphs for <f> and A y with J/T c 2 = 0.0094,/ioo = 5, 
L = 4, e = 0.6, and a = 0.5, are shown in Fig. 1. 

Then, in Fig. 2, we show that the DC current of p-wave junction is proportional to the 
sine of the phase difference across the junction, that is to say, the red dots coming from the 
numerical calculations match with the the black solid sine curve. Fitting the sine curve to 
the data, we can obtain the maximum current across the junction: 

Jmax/Tf « 1.078. (3.1) 
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Figure 2. Supcrfluid current Jmnx/T^ as the function of the phase difference 7. The black line is 
the sine curve. The parameters are set to [x^ = 5, L = 3, e = 0.6, = 0.5. 

The dependence of J max on the width of the gap is shown in Fig. 3(a). The graph 
predicts an exponential decay with the growing width of the gap in J max : 

^ = A e~i, (3.2) 

where £ is the coherence length. In Fig. 3(b), the condensate (O) = ip^ at zero current is 
shown. The graph also predicts an exponential decay with the growing width of the gap 
in t/ m ax' 

T2 \j=o = A 1 e 2 «. (3.3) 

Fitting Eq. (3.2) and Eq. (3.3) with the two sets of data, we can obtain {^,Aq} ps 
{1.11,16.42} and {^,^1} fa {1.15,43.89} for Eq. (3.2) and Eq. (3.3), respectively. The 
disagreement of two values of £ is about 4 percent. In Fig. 4, we obtain the relation of 
J max and T. One can verify that near the critical temperature T c , J max can reach zero. 
Since e = 0.6, we show the region corresponds to T/T c < 0.6, which depicts the character 
of a p-wave Josephson junction. 

4 Conclusion 

In this paper, we construct a holographic model for p-wave SNS Josephson junction with 
DC current. In the background of (3+l)-dimensional gravity, we solve a set of couple, 
partial differential equations of a SU(2) gauge field numerically. By choosing spatial y- 
dependence /j,, we find the the DC current of p-wave junction is proportional to the sine of 
the phase difference across the junction. Moreover, the graph which predicts an exponential 
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Figure 3. The curve of J max and (O) x= o on L. The parameters are set to /i^ = 5, e = 0.6, 
a = 0.5. 

decay with the growing width of gap in J max is obtained. At last, the curve of J max and T 
is also obtained. One can see that the model of the holographic p-wave Josephson junction 
in our work can match precisely with the results for condensed matter physics. In future, 
It would be interested for us to extend to investigate the holographic d-wave Josephson 
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Figure 4. The curve of J max on T, The parameters are set to e = 0.6, a = 0.5 



junction. 
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